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1 Introduction 


If the set of minima At of the potential energy of a field is not 2-connected, for example 
if it is a 2-sphere S'^ or real projective space MP^, then it admits solutions called global 
monopoles. In general the total energy of a global monopole exhibits a long distance 
divergence (linear divergence). 

In the literature the global monopole has been considered mostly in a cosmological 
context where the gravitational energy, their Goldstone boson radiation and gravitational- 
wave emission and finally pair annihilation sparked the primary interest dlEllal. Global 
monopoles may also be realized in nematic fluids jlj. 

The stability of a single global monopole has been debated in the early literature. 
It was argued that the global monopole could itself collapse to a zero-energy solution 
[S], but it was shortly after argued not to hold true [H]. Ref. [7j further argued that the 
observation of Ref. [5] should not be interpreted as a collapse but rather as an acceleration 
of the entire monopole. The false conclusion was thus reached by artificially holding the 
origin of the monopole fixed. The bottom line is that the single global monopole is indeed 
stable, both with respect to radial and angular perturbations [7]. 

Ref. [7] made a further study of the mutual forces between two different types of a 
global monopole and an anti-monopole, which we will call type B and type C. We find - 
in agreement with Ref. [7j - that the type C monopole and anti-monopole attract each 
other and the type B repel each other. We find, however, different values of the forces 
than Ref. [7j does. We further generalize this calculation to two monopole-monopole cases 
which we will call type A and type B and they are both mutually repulsive. 

The literature has nothing much to report about global multimonopole configurations. 
This is because they are unstable. For completeness we carry out a numerical calculation 
of a perturbed two-monopole configuration in [Appendix A| to confirm that they are indeed 
unstable; that is they repel. 

We contemplate an interaction term which may be added to the global monopole 
Lagrangian in order to stabilize said system. This is argued, however, to be problematic 
due to the possibility that the monopoles just to shrink upon such an intervention. We 
have also checked this with numerical calculations. 

As a different approach on the path to stabilize the global two-monopole configura- 
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tion, we consider a modification of the global monopole Lagrangian more similar to a 
nontrivial effective metric. We choose to construct said effective metric being a func¬ 
tion of a Skyrmion field. The stabilization is thus dynamic but not renormalizable. The 
Skyrme Lagrangian itself is not renormalizable either, so we do not regard that as 

a big problem but the model as an effective field theory model generated by some fun¬ 
damental theory. A single gauged monopole inside a U(l) gauged single Skyrmion has 
been studied in the literature already ng. We are studying the global analog of such 
a configuration and the multimonopole instead of the single monopole. We show that 
the global two-monopole can indeed be stabilized inside a Skyrmion for certain choices of 
parameters. 

The paper is organized as follows. In the next section we review a single global 
monopole to remind the reader and set the notation. In Sec. we review the monopole- 
anti-monopoles interactions in the point-charge approximation and make the same type of 
analysis for the two-monopole system. After Sec. we no longer use the point-charge ap¬ 
proximation and all the calculations and considerations are carried out with the monopole 
profile functions and hence with true equations of motion. In Sec. we discuss the pos¬ 
sibility of stabilizing the global two-monopole using an interaction term. In Sec. we 
construct the two-monopole inside a Skyrmion and dub it the Skyrmopole. Finally, we 
conclude in Sec. The [Appendix A| presents a numerical calculation showing the insta¬ 
bility of the normal global two-monopole. 

2 A single global monopole 

The theory consists of an adjoint-valued scalar field, d> = the gauge indices 

a = 1,2,3 are summed over, cr“ are the Pauli matrices and the Lagrangian density is 

£ = -^(l-TV|*^|)t (2.1) 

which has the equation of motion 

= -\v^ (1 - Tr[d>2]) $. (2.2) 

The Greek letters are used for spacetime indices, fi = 0,1,2,3 and we use the mostly- 
positive metric. The potential breaks the global SO(3)-symmetry down to SO(2) and the 
coset SO(3)/SO(2) ~ when mapped from spatial infinity, which is also S^, defines the 
topological charge by means of the second homotopy group. The topological charge of 
the monopole is given by the winding number, i.e. 

Q = -^f (2.3) 
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where a,b,c = 1, 2,3 are gauge indices and i,j = 1, 2,3 are spatial indices. The charge Q 
is also called the monopole number. The mass of the monopole is the total energy in the 
rest frame 



(2.4) 


Let us hrst review the standard hedgehog conhguration which describes a single global 
monopole in the theory at hand. Since Q = 1, we can use a spherical Ansatz 


$ = -/j(r)x“cr“, 
r 


for which the energy density reads 


V = - It, (2.6) 

2 rpZ 

where ' = dr- Defining a dimensionless coordinate, p = \/Aur, we can write the energy 
density as 

= + + (2.7) 


where h = h(p); now ' = dp and the equation of motion reads 

h"+-h'-l)h = {). (2.8) 

P P 

The boundary conditions for the monopole are: h(0) = 0 and h(oo) = 1. Expanding 
around p = 0, we get 

h = ap -—p H p +0{p), (2.9) 

where a > 0 is a constant (also called the shooting parameter). This tells us that the 
energy density around the center goes as 




( 2 . 10 ) 


i.e. roughly constant for small values of a. Asymptotically, h is very close to one, so 
linearizing about /? = 1 — h, we get the asymptotic equation of motion 


/5" +V-4(/5-1)-2/5 = 0, 

P P 


which has the exact solution 


h = 1 — B = 1 


( 2 . 11 ) 


( 2 . 12 ) 
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This in turn determines the asymptotic energy density as 


11 1 2 1 


(2J3) 


which for large p goes like 1/p^ and hence the total energy diverges linearly. Finally, the 
topological charge can easily be calculated using Eq. (2.3) as 


Q = 

^ Stt 


dS*’ = 1. 


(2.14) 




(a) (b) 



(c) 


Figure 1: (a) profile function h, (b) energy density and (c) energy density on loglog scale 
for the hedgehog ((5 = 1) monopole. 


In Fig.[^is shown the profile function and energy density of the hedgehog configuration. 
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3 Charge-two and charge-zero monopoles in the point-charge 
approximation 


Let us now turn to the double-winding global monopole configuration. The overall 
orientation of the monopole is not important, but the relative orientation is important 
for the multi-monopole solution. In the following we will write the scalar field as 

$ = h {sin 5 -cos fc, sin p sin A:, cos p} , (3.1) 


which describes a charge-two monopole if k winds twice and g winds once. Another 
possibility is to make g wind two times and k one time, but it requires a twist (sign 
change) of for instance k, in order not to create a monopole-anti-monopole pair. To see 
the charge explicitly, let us evaluate the integral (2.3) for g = m6 and k = n<p oX some 
asymptotic distance where = 1: 


Q = — 

^ Stt 


= nsin^ — 


m7T\ 


- 


2 / 


(3.2) 


We will study a number of cases in turn. 

For physical and everywhere regular global monopoles, the profile function h cannot 
take the value 1 over all space, but should have at least Q zeros, counted with multiplicity. 
Regularity dictates that the zeros of h coincide with the origins of where the phases of g 
and k wind about. The approximation of setting h := I can physically be seen as the limit 
of very small monopole sizes compared to the separation distances. Once we set h := 1, 
we introduce singularities at these origins or positions of the charges. These singularities 
are dealt with by subtracting the same singularity of a single singular global monopole 
field. The difference is then a regular field. This approximation captures the asymptotic 
interaction of global monopoles, where the potential has no effect, but only the winding 
phases interact at long distances via the kinetic term. This interaction is important as it 
can never be neglected. 

Independent of the local regularity of the global monopole, the total energy of the 
configuration always diverges linearly. This divergence is just an artifact of putting a 
global monopole in an infinite space. In real world applications, such as liquid crystals, 
global monopoles only live in a finite region of space. In this paper we will consider the 
global monopoles on an infinite space, although the numerical calculations will of course 
only be carried out on a finite space, just like the physical global monopoles enjoy. 

A comment is in order about the Ansatze that we will use for the function g and 
k in the next subsections. They are not describing solutions, but describe the winding 
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that necessarily will be present on a given timeslice in a would-be solution. Logically, 
there are two possible situations. If the two charges attract, a solution will be close to 
a configuration where the charges are coincident. By making an educated guess for the 
function /j, the full, time-dependent solution can then be found using the full equations 
of motion. The other possibility is that the charges repel each other. In that case, no 
stable solution exists and this type of configuration is never a static solution, but at 
best a snapshot of a time dependent field configuration that evolves in time. Real time 
dependent solutions to the equations of motion can be constructed, but since we are not 
interested in unstable time dependent configurations, we will simply disregard such cases 
as being what they are: unstable. 


3.1 Type A 


The first case is the charge-two monopole where k winds twice and g once 
^fuii = h {sin 6 cos(0i ± 02 ), sin 0 sin(0i ± 02 ), cos 9} , 


where we define the angles as 


6 = arccos -, 0 i = arctan —-, 
r X — a 


02 = arctan 


y 


X d 


(3.3) 


(3.4) 


and the two monopole charges are separated by a distance 2d. We will treat a case of a 
monopole-anti-monopole on the same footing as the monopole-monopole by choosing the 
lower sign. The monopole charge of the configuration does not depend on d; hence if we 


set d = 0, we can see from Eq. (3.2) that this configuration has charge 2 for the upper 


sign and charge zero for the lower sign. 

In order to study the stability of this monopole, let us make a point-charge approx¬ 
imation, where h := 1. This approximation corresponds to ignoring the potential and 
taking into consideration the asymptotic effect of two charges affected only by the kinetic 
term. See also the discussion in the previous subsection. 


The form (3.3) has been chosen for its simplicity. However the price of this simple 
form is that, when d 7 ^ 0 , the field •hfuu does not wind around the vacuum manifold when 
restricted to a small sphere centered on the singular point at x = ±d. Indeed, on such a 
small sphere the value of cos 9 is always near zero. The solution restricted to a 2-surface 
can only have nontrivial winding on the vacuum manifold if at some point sind = 0 , 
which means that the 2 -surface must intersect the axis and so extend a distance d from 
the singular point. In fact, the kinetic energy diverges along the infinite lines extending 
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along the axis from the two singular points, which leads to a short-distance logarithmic 
divergence in the energy. 


This divergence is simply an artifact of our simple functional form (3.3). In Appendix 


[A| and in Sec. when we consider true solutions of the equations of motion, the field 
configurations will be regular everywhere and so will the energy be. However, in the 
present section we are only interested in the interactions between monopoles. What we 
will do is to subtract off a counter term, which cancels the divergence. We will call this 
difference the interaction energy and it is free from short-distance divergences. A concern 
may be that the exact form of the counter term may change the functional behavior of 
the interaction energy; however, to leading order in d it is independent of the counter 
term. To see this, we define the following interaction potential 


= £^[$full] - ^[$l] - E[^2\, 

where is the z-th single monopole which we subtract off 

d>j = h {sin 6 cos 0*, sin 9 sin 0,, cos 6} . 


(3.5) 


(3.6) 


The single monopoles that are subtracted have the same two fines of divergences at a: = 
±d, y = 0 and so the respective contributions of their short-distance divergences to the 
interaction energy cancel. The interaction energy is sufficient to determine the radial force 
between two monopoles, which is simply its derivative with respect to d. 
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The interaction energy is 


E. 


int 


d^x 


- 2d‘^ [x^ - y'^ T + a;^)] + (1 T 2)(x^ + y'^f 

[{d — xy + 1/2] [(d + x)2 + //2j (^2^2 _|_ y2 _|_ 


J p<L,\z\<L 

1 , d‘^ ± 2 (i^p 2 + (1 ip 2)p^ - 2 d?p^ cos 20 L 

dp I d(p -— nj 9 -arctan I — 


'0 JQ 
cL 


— 2(^2 p2 cos 20 


= —27r J dp ^ ± sign((i — p) 


2p2 


d? + p^ 


arctan 


P 


1T 2 , , ^ L 2 7 2c?L 

--—7r(7r + 2 log 2)L Sttci arctan — ± 27r a arctan —-—— 

z cl L/^ — 

f d^\ L -\- d 

=F47rLlog ( 1 + ^ ) T47rdlog21og ^ ^ 


± -vrdLio 
2 


d-L 
d + L 


=F ^vrc?Li2 


d + L 

d-L 


^ t / t - ^ + ^ 
± 7rdLi2--7 

d — L 


d — L ISO! ■ {d + L)(d + iL) i<^t ■ {d — L)(d + iL) 

T 7rdLi2 -7 ± 7r(i3?Li2 y-. -irfy-- ttt T 7rd3fJLi2 - 


‘d + L 


{d-L){d-iL) 


{d + L){d-iL) 


■ ^ ^ ^7 r(7r + 2 log 2)L 47r^d ^ ) • 


(3.7) 


The domain of integration here is chosen, for convenience, to be a cylinder (as opposed 
to a sphere). The first term in the linear expansion does depend on the geometry of the 
integration region (cylinder versus sphere etc.), however the second term, corresponding 
to the force does not. Again the upper sign is for the monopole-monopole case and the 
lower sign is for the monopole-anti-monopole case. 

The mutual force between the monopole and the (anti-)monopole for small d is 

\Fmt ±47r^ (3.8) 

and so the two monopoles (monopole-anti-monopole) repel (attract). 

Fig.g shows the interaction energies for both cases of type A. 


3.2 Type B 


The next case is a charge-two monopole where k only winds once and g winds twice [I] 
^fuii = h {sin(di 62 ) cos 0, e sin(di -h 62 ) sin cf), cos(di + ^ 2 )} , (3.9) 





















(a) monopole-monopole 


(b) monopole-anti-monopole 


Figure 2: Interaction energy of (a) the type A monopole-monopole pair and (b) the type 
A monopole-anti-monopole pair as a function of separation distance 2 d for a cut-off value 
L = 20. The blue dashed line shows linear expansions in d. 


where we define the angles as 


9i = arccos 


z — d 


+ {z — dy^ 


z -\- d y 

Uo = arccos — , =, ©i = arctan —, 

^/x^+y^ + (z + d)^ x' 

(3.10) 


and the monopole charges are separated again by a distance 2 d. e is the twist function; if 
e = 1 the configuration describes a monopole-anti-monopole pair, which was considered 
in Ref. [7]. In order to get a charge-two monopole configuration, we need to twist the 
second cycle in q, i.e., 

e = (-l)LfJ, (3.11) 

where [x\ is the floor function of x and 


6 = arccos 


^yx'^ + y'^ + z'^^ 


(3.12) 


is 6 measured from the origin which is the midpoint between the two monopoles. 


Let us now study the point-charge approximation, where h := 1. 
interaction potential (3.5) and $j, which is the f-th monopole, is now 


We again use the 
given by 


= h {sin 6 i cos 0, sin 9i sin 0, cos dj} . 


(3.13) 
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The interaction energy thus reads 


9 -^int 


(fx 


(f 


r<L ~ 


^ dz 
2 d J_i^ 


(f-z 


2 r (d-zf , d^ + L‘^ + 2dz 

log 7TT-X^ + log 


ttL / 27r^d 


(d + 2 )^ 


ttL fL^ 2L 3d\ fd 

+ t U + t-i: I 


^ , d\ / 2dL 


8 ttL n {d^ 

—-*» + o - 


(3.14) 


where 

?7(x) = Li2(x) - Li2(-x), (3.15) 

and Li 2 is the Spence function (also called dilogarithm). Note that the negative sign in 
front of the order d term in the expansion means that the system has negative binding 
energy and the mutual force is repulsive 

i-^int-TT^- (3.16) 


Note that the interaction energy remains the same for both the case with the twist 
function (the monopole-monopole case) and without the twist function (the monopole- 
anti-monopole case). 

In Fig. is shown the interaction energy of the charge-two monopole (which is the 
same as the monopole-anti-monopole) of type B. 

Let us mention that this configuration has no lines of divergences, so only the singular 
points at the origins of the charges are still present in $fuii but Ei-^t is again free from 
these singularities. 

We will now consider the total energy, which is the same for both the monopole- 


10 












Figure 3: Interaction energy of the type B monopole-monopole pair (which is the same as 
the type B monopole-anti-monopole) as a function of separation distance 2d for a cut-off 
value L = 20. The blue dashed line shows the linear expansion in d. 


monopole and monopole-anti-monopole case, 

-^total 


_|_ ^2 


= 2 / d^T — 

Jr<L {<P + 


- r ^ 

J-L z 


— 2 dz) log 


{d - zf 


ttL f 2'K‘^d 


d^ + L^- 2 dz 
ttL lOL lid 


-h — z^ + 2 dz) log 


(d -|- z') 


d2 + L2 + 2 dz 




d 


. arctanh ( y 
J-J / \ J-J 


2'Rdri ( — 


Tvdrj 


2dL 

d^Tl^ 


dTT" + C» ( ^ 


_ 32wL 

~ 3 ' ■ V L 

Again the total energy £'[d>fuii] diverges linearly in L for both configurations. 


(3.17) 


3.3 Type C 

The last case, which we consider for completeness, is a monopole-anti-monopole where k 
winds once and g winds and unwinds at spatially separate locations [7j 


^fuii = h {sin(di — 62 ) cos 0, sin(di — 62 ) sin (p, cos(di — ^ 2 )} , 


(3.18) 


where the angles are given in Eq. (3.10) and the monopole and anti-monopole charges 


are separated again by a distance 2d. This configuration was also considered in Ref. [7j. 
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We again use the point-charge approximation where h := 1. The interaction potential 
is given in Eq. (3.5) and the single monopole fields are the same as in Eq. (3.13). The 


interaction energy thus is 
d^x 


— 2 


d? -r‘^ 


r<L {(P + 

dz 


= TT 


l-L 


Z 


, , {z — df {z + dY 

{z - d) log + {z + d) log 


d? + L'^ — 2dz 


d? + + 2dz 


'KiLd — id + Vd . j. A 1 f d . 

- d - - \d —L)2 ~ ~ \ L ~ 

d? 


2dL 
d^ + L^ 


= —SttL + 2d7r'^ + O \ — 


(3.19) 


Note that the positive sign in front of the order d term in the expansion means that the 
system has a positive binding energy and the mutual force is attractive 


^ -27r'^ 


(3.20) 


In Eig. I^is shown the interaction energy of the monopole-anti-monopole of type C. 



Eigure 4: Interaction energy of the type C monopole-anti-monopole as a function of 
separation distance 2d for a cut-off value L = 20. The blue dashed line shows the linear 
expansion in d. 
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Finally, we will evaluate the total energy for the type C monopole-anti-monopole 

— £;total = 4^2 d?x ^^2 + ^ 2)2 _ 4 ^ 2^2 

■ id^zf , d? + L^-2dz 
{d-zf ^ d? + L^ + 2dz_ 

= 2d7r^ + O . (3.21) 

Notice that this and only this case (of those considered here) has a hnite total energy. 

Let us hnally mention that the interaction energy is regular in this case but has a 
linear divergence in the integral (the constant), whereas the total energy has 2 singular 
points at the monopole charges but has a finite total energy. 



3.4 Summary 

The conclusion we can draw from the point-charge approximation is that the charge-two 
monopole is unstable given both orientations we have constructed. The mutual repulsion 
could be slightly altered by including the profile function h of the configuration, which in 
turn includes the effect of the potential. Asymptotically, however, the instability of the 
global charge-two monopole remains, as we have demonstrated above. For concreteness 
we carry out a numerical calculation of a perturbed two-monopole configuration in [Ap¬ 
pendix A| to confirm that they are indeed unstable. In the next section we will begin to 
contemplate how to stabilize the global two-monopole. 


4 Caging global monopoles 

In this section and in the remainder of this paper, we consider the monopole profile 
function and hence the full equations of motion. The first attempt at stabilizing the global 
monopoles is to modify the potential in such a way that it is energetically favorable to 
remain near the spatial origin of a certain coordinate frame. This origin could then be 
determined dynamically by a solution to a different sector, which we here for concreteness 
take to be a one-Skyrmion sector. 
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In this section, we wonld like to advocate that it is not as trivial to realize this idea 
as one might expect. We consider for concreteness the following theory 

A?d 9 

-£ = -TV(a„*)= + —(1 + "4)“ (1 - rvM) , (4.i) 

where the field n 4 is a component of the Skyrmion sector which obeys the boundary 
conditions 714 ( 0 ) = — 1 (at the origin of the coordinate frame) and 774 ( 00 ) —>■ 1. 

The idea is that at the origin, the potential vanishes and hence the contribution to the 
energy from the potential is minimal when the monopole is situated at the origin. If the 
monopole tried to get out, it should start to feel the potential growing and thus it would 
be energetically favorable to return to the origin. The hope is then that by tuning the 
potential parameters, the two monopoles still prefer to stay at the origin although they 
are mutually repulsive. There is however a monkey wrench around the corner. Hence, 
instead of presenting a numerical calculation, let us examine the expectations with a very 
crude calculation. 

Let us assume that the Skyrmion has the size Lgk, which for all practical purposes 
here just means that 

2?^ 

774(r) = -1 + —, for r < L^k, (4.2) 

r^sk 

and r is the radial coordinate measured from the origin. Let us further assume that a 
single monopole has the profile shape = a^'r/L^ for r < L^. We consider the case 
Lsk ^ Lm, which is what we would prefer for physical reasons. The repulsive force comes 
from the kinetic term and was estimated in Sec. |3]to be 


-^repulsive 477 V , 

and the attractive force due to the above potential can be estimated as 

^2l^2 


(4.3) 


\.,4 

-^attractive = j ^[1 + n^{d)T (l “ Tr[$^]) 

2q;A?;^ 


L 


dr I 


sk 


O'—! / 2 \ 2 

fi- — 

Lsk J \ L^ 


(4.4) 


Since we assumed that the monopole size is much smaller than the size of the Skyrmion 
(Z/m Lsk), we can crudely pull the factor of (2d/Lsk)“~^ out of the integral and we get 


IGctTrAn^L^ ( ‘^d\ 

-^attractive TrTrT I ~r I 

105Lsk V-^sk/ 


a—1 


(4.5) 
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Now, the distance from the center of the origin d, has to be small for the potential to be 
effective (in the sense of exerting a force on the monopole); that is, d < Lskj^ Let us now 
compare the forces, i.e. the attractive one versus the repulsive one 


AaXv^Ll^ 



a—1 


lOSZ/gi^ \Lq\^ 


> TT. 


(4.6) 


The inequality shows that a larger Skyrmion provides a smaller attractive force (this is 


also expected as the derivative of the field profile becomes smaller). The real kicker is that 
by trying to increase the potential factor, i.e. Xv^, the monopole size decreases by the same 
factor, namely ~ Xv^. This means that the only parameter we can dial to stabilize 
the monopole(s) is a. However, a large value of a also makes the factor (2d/Lsk)"~^ small. 
In other words, once we try to squeeze in the monopole(s) by the potential, they shrink 


and can escape as easily as before. 

This is not supposed to be a proof of no-go, but merely a justification of not construct¬ 
ing the theory with just a modified potential. We have tried several numerical attempts 
and have observed exactly this behavior, namely that the monopoles shrink and are able 
to escape. We do not consider those attempts as a proof of no-go either and will not 
present this study here. In the next section we will pursue a different strategy. 

5 Skyrmopole 

In this section we will place a global two-monopole inside a potential so that each 
monopole will not run away despite mutual repulsion. In view of the estimates of the 
last section, we consider a coupling between a one-Skyrmion and the two-monopole via a 
metric-like prefactor of the monopole part of the Lagrangian density. The aim is that the 
monopole constituents prefer to stay in the center of the “potential”, so much so that the 
repulsion is negligible. If a monopole constituent is to exit the Skyrmion, it will enjoy the 
instabilities described in the Sec. IH 

Let us first contemplate the construction in order to know what to expect. We start 
with a Skyrmion which has a kinetic term, a Skyrme term and a mass term. The size 
of the Skyrmion is given by the balance of these three terms. The global monopole, on 
the other hand, has a size squared given by the inverse coefficient of the self-interacting 

^One may naively think that the factor is helping at making the force bigger when d > Z/sk/2, but 
in the real Skyrmion solution, the tail is much flatter at such distances from the origin and hence the 
derivative becomes small. 
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potential term. In front of the monopole sector we multiply by a factor depending on the 
profile of the Skyrmion sector. 

Putting the pieces together gives us the Lagrangian density 

= -^sk + (5-1) 


where the Skyrmion Lagrangian is simply 


-jC-sk = ^(9^^ • ^ [(a^n • • d^Tof] + m^(l - 724), 

the monopole Lagrangian reads 

- £„ = ^ (1 - Tr[$ 2 ])' , 

and the coupling is chosen to be 



( 6 . 2 ) 


(5.3) 


(5.4) 


where the fields Up, with p = 1 ,... ,4 are unitless 0(4) fields describing the Skyrmion 
part of the configuration and •fi = 'I>“( 7 “, a = 1,2,3 are SU( 2 ) adjoint unitless monopole 
fields. The coefiicients in the Lagrangians are C 2 and m with units of mass squared, v has 
units of mass, and finally C 4 , A, b and a are unitless parameters. C 2 , C 4 and m determine 
the size of the Skyrmion and y/Xv determines the size of the monopole (s). b is included 
in order for the monopole energy not to vanish at the center of the Skyrmion and larger 
values of a translate to a larger attractive force among the monopoles. One may interpret 
the coupling, Q, as an effective metric felt by the monopole sector of the theory (although 
here it is just incarnated as a non-renormalizable field theory coupling). 


If we compare to the considerations made in the last section, the difference is that now 
the size of the monopoles and the size of the energy barrier are independent parameters 
of the theory. 


The static equations of motion are 


C2d‘^np + C4 (9jn • 9*11) d‘^np + C4 {didjn ■ dju) diUp — C4 (cl^n • din) diUp 


av"^ ^ 1 + 6 + 724 Y 

o^—l 


+ 772^6^^ = 0, 

(5.5) 

2+6\ 2+6 / 



ry 

+ , 9.7249'^$“ + At;^ (1 

1 + 6 + 77.4 

Tr[d>2]) $ = 0. 

(5.6) 

The vacuum equations are 


Tr)*!*^] = 1, 724 = 1. 


(5.7) 
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Note that even without the mass term for n, there is no explicit breaking of 0(4) symme¬ 
try, the symmetry is still broken spontaneously down to 0(3) for any finite energy state. 
This breaking, 0(4) —?• 0(3), supports the Skyrmion(s). The spontaneous breaking of 
SU(2) down to U(l) by the potential supports the monopoles. 


There is a subtlety due to the coupling of the monopole theory to the Skyrmion sector 
in this way, namely that the second last term in Eq. (5.5) cannot exceed the value of the 
mass term in order for the field, n 4 , to be able to return to its vacuum expectation value. 
When this is not fulfilled, so-called pion condensation occurs and the vacuum equations 
are completely different, see e.g. Ref. m- This leads to the following (strong) condition: 


av"^ 
2 + b 


lTr(8„tf + ^ (1 - Tr|J.2])" 


< m^. 


(5.8) 


where we have neglected the factor of ((1 -|- 6 -|- n 4)/(2 -|- 6))““^, as it is always smaller or 
equal to unity (by construction). This means that the region in parameter space we are 
interested in is where a is big (to get a sizable attractive force) and v is small (to avoid 
breaking the vacuum conditions)]^ 

The Skyrmion charge is also a topological charge and can be calculated as 

B = J e^'^^'^'e^^’^dinpdjnqdknrns. (5.9) 


In order to get a feeling for the parameters in the theory, let us consider a 1-Skyrmion 
with a single monopole inside, i.e. Q = B = 1. For this we can choose the following 
Ansatze 


n= sin/(r), sin/(r), sin/(r), cos/(r)} , (5.10) 

d> = h(r)a;“cr“, (5-11) 

where £“ = x“/r is the spatial unit vector and r = \x\ is the radial coordinate. Inserting 
the above into the Lagrangians (5.2) and (5.3), we get 

-Ak = ^fr + ^sinV + ^sin2(/)/r + ^ sin^/ +772^(1 - cos/), (5.12) 


-Q{nA)C.^ = 


1 -|- 6 -|- cos / 

2Tb 




(6.13) 


^If the boundary conditions are switched so that 77.4 = 1 at the origin and 714 = —1 at spatial infinity, 
the same condition holds with 77.4 ^ —77,4, i.e. it does not ameliorate the problem. 
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which give rise to the following equations of motion 


C 2 ( frr H- fr -^ siu 2/ ) + C4 


4 sm^{f)frr + 4 sin(2/)/^^ - ^ sin^ / sin 2/ 


— sin / + 


an 


2 + 6 


1 + 6 + cos / 

2 T 6 


Q:—1 r -1 -| \ 2 

hl + -h^ + ^il-h^f sin/ = 0, (5.14) 
2 4 


J ^hr + An2(l - = 0, (5.15) 

r r"' 1 + 6 + cos / 


where = drf. Let us now expand the equations of motion around r = 0 as 


/ = TT + /ir + 


h = hir + + 


(5.16) 


The equations of motion to leading order in r determine /2 = ha = 0. We can thus 
calculate the energy density around the origin as 


^ 3 ,2 3 ,4 o ^ 2 f b 

£ = y2/. + 2C4/,+2m + 2 +^ 


A? + hi* 

o 


+ 0 +). 


(6.17) 


Using instead a scaling argument due to Derrick lEI, we can estimate the size of the 
Skyrmion as follows. Neglecting the monopole fields (as we assume that -C (c 2 ,m)), 
we denote by eaq the kinetic energy terms with 2 and 4 derivatives and u is the unitless 
potential 

U = caCa + €464 + (5.18) 

Rescaling now the spatial coordinates x* [ix\ we get 


r. 1 la 

E -caca + ^^0464 H —-m u, 

[I /j.'^ 


(6.19) 


which has an equilibrium when 

/i^C 4 e 4 — ^ 1^0262 — 3m^u = 0. (5.20) 

If we neglect the kinetic term, assuming that the mass is big and the size of the Skyrmion 
is also big, then the characteristic length scale and hence the Skyrmion size is roughly 

(5.21) 


L. 


sk 


C 4 

3m? 


The hnal step is to identify fi in the expansion of the solution as an order-one number 
times the inverse length scale: fi ~ 1/Lsk- Now we can write the energy density of the 
Skyrmion part as 

3V3c2m 2 
^sk ^ -^ + 5m , 


2 ^/^ 

whereas the energy density of the monopole part is roughly 


^rr 


l^'’\2 + b 


(5.22) 


(5.23) 
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5.1 A single Skyrmopole 


Let us try to solve the coupled equations numerically. We use the finite-difference method 
on a cubic lattice of sizes 121 ^ and 201 ^ and evolve a linear time operator, which is standard 
in the relaxation method]^ After the equations of motion are satisffed locally better than 
the one-permille level we stop the relaxation algorithm. We use the Ansatze (5.10)-(5.11) 
with appropriate radial profile functions as the initial configuration for the relaxation 
algorithm. We have chosen the parameters as C 2 = 1, C 4 = 4, m = 1, u = 1/4, A = 128, 
6 = 1/2 and a = 4 in order to satisfy the constraint ( |5.8[ ). In Fig.|^is shown a single global 
monopole {Q = 1) inside a single Skyrmion (B = 1). The figure shows the 3-dimensional 
isosurfaces of the Skyrmion charge and monopole charge densities, respectively, each at 
their respective half-maximum values. 


The size and density (lattice spacing) is the same as that used for the two-monopole 
inside a Skyrmion and we have optimized said dimensions for the two-monopole. This 
is why we do not capture the Q = 1 monopole charge better than to about three and a 
half percent on this lattice. As we will see shortly, the two-monopole charge is far better 
contained. 


In Fig. [^are shown xy slices at z = 0 of the monopole charge and Skyrmion charge, 
respectively. Both have the shape of solid balls and the monopole has roughly half the 
diameter of that of the Skyrmion. 

In Fig. are shown six xy slices at z = 0 of energy densities. More precisely of the 
monopole kinetic energy, the Skyrmion kinetic energy, the monopole potential, the pion 
mass term (Skyrmion potential), the total energy and finally the logarithm of the total 
energy. The shape of the Skyrmion is ball-like while that of the monopole is rather like a 
hollow sphere. This is in contradistinction to the monopole charge-density distribution. 
The reason is that the monopole energy is multiplied by a factor Q, which amplifies the 
energy outside of where the monopole charge is located. This is no accident. This is made 
by design of the function Q in order to provide an attractive force on the two monopole 
constituents. One can think of this construction as making an energy barrier that the 
monopoles do not want to cross and hence choose to stick with each other instead (even 
though they are repellent to one another). 

The reason why this calculation is so numerically challenging is that the monopole 

^Clearly, the single Skyrmopole case, Q = B = 1 , can be solved using just the ordinary differential 
equations and not the three-dimensional partial differential equations (PDEs). In fact we solved this case 
also using an ODE solver and compared the solution to the one shown here found by solving the PDE as 
a check on our code. 
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Figure 5: A global monopole inside a Skyrmion. The figure shows isosurfaces of the 
monopole charge and Skyrmion charge densities, respectively, at half their respective 
maximum values. The gray cloud is the Skyrmion charge isosurface and the colored 
isosurface is the monopole charge. The coloring is made using an HSL (hue-saturation- 
lightness) map from the monopole field such that arg($i/'F 2 ) is mapped to the hue and 
I$ 3 1 determines the lightness. The parameters are chosen as: C 2 = 1, C 4 = 4, m = 1, 
n = 1/4, A = 128, b = 1/2 and a = 4, and the numerically integrated charges are 
B = 0.99993 and Q = 0.96707. The calculation is made on a 201^ cubic lattice. 
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(a) (b) 


Figure 6 : Topological charges for the global monopole inside a Skyrmion; (a) is the 
monopole charge and (b) is the Skyrmion charge. The figures show slices in the xy plane 
for z = 0. The numerically integrated charges are B = 0.99993 and Q = 0.96707. 


charge is much smaller than the corresponding energy distribution. This is because the 
prefactor makes the majority of the monopole energy density sit outside of where the 
monopole charge is situated. This separation of scales requires both very large length 
scales (for capturing the energy) and a very dense lattice for capturing the charge at the 
center. 


5.2 The type A, Q — 2 Skyrmopole 


We are now ready to compute the two-monopole inside a single Skyrmion; we can call it 
the Q = 2 Skyrmopole. As the initial condition we use the Ansatze (5.10) and (3.3) with 
d = 0. We take the profile function of the initial condition to be 


h = tanh Ikt] 


where r is the radial coordinate and k is an appropriately chosen constant. We use the 
same lattice and same parameters as in the previous subsection, i.e. C 2 = 1, C 4 = 4, 
m = l,n = l/4, A = 128, 6 = 1/2 and a = 4. We again satisfy the constraint ( |5.8| ). In 
Fig. I is shown the global two-monopole [Q = 2 ) inside a single Skyrmion (5 = 1 ). The 
figure shows the 3 -dimensional isosurfaces of the Skyrmion charge and monopole charge 
densities, respectively, each at their respective half-maximum values. 
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(d) (e) (f) 


Figure 7: Energies for the global monopole inside a Skyrmion; (a) is the kinetic term 
of the monopole (with the prefactor Q), (b) is the Skyrmion kinetic energy, (c) is the 
potential for the monopole (again with the prefactor Q), (d) is the pion-mass term of the 
Skyrmion sector, (e) is the total energy of the configuration and (f) is the logarithm of 
the total energy. The figures show slices in the xy plane for z = 0. 
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Figure 8 : The global two-monopole of type A inside a single Skyrmion. The figure shows 
isosurfaces of the monopole charge and Skyrmion charge densities, respectively, at half 
their respective maximum values. The gray cloud is the Skyrmion charge isosurface and 
the colored isosurface of the shape of a torus is the monopole charge. The coloring is 
made using an HSL (hue-saturation-lightness) map from the monopole field such that 
arg($i/'F 2 ) is mapped to the hue and |$ 3 | determines the lightness. The parameters are 
chosen as: C 2 = 1 , C 4 = 4, m = 1 , u = 1/4, A = 128, b = 1/2 and a = 4, and the 
numerically integrated charges are B = 0.99993 and Q = 1.9863. The calculation is made 
on a 201 ^ cubic lattice. 
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Figure 9: Topological charges for the global two-monopole inside a single Skyrmion; (a) 
is the monopole charge and (b) is the Skyrmion charge. The hgures show slices in the xy 
plane for z = 0. The numerically integrated charges are B = 0.99993 and Q = 1.9863. 


In Fig. [^are shown xy slices at z = 0 of the monopole charge and Skyrmion charge, 
respectively. The Skyrmion is still a ball while the charge distribution of the two-monopole 
has the shape of a torus. Notice that the two charge distributions have almost the same 
size. 


In Fig. are shown six xy slices at z = 0 of energy densities. Again the subfigures 
show the monopole kinetic energy, the Skyrmion kinetic energy, the monopole potential, 
the pion mass term (Skyrmion potential), the total energy and finally the logarithm of 
the total energy. As was the case of the shapes in the one-monopole case, the Skyrmion 
energy is ball-like while the monopole energy has the shape of a hollow sphere. 


We numerically integrate the energy densities shown in Fig. [T^ and give the results 
in Tab. [!} We integrate only the solution inside the box on which we find the solution, 
which has the size 20^. This captures the total energy of the Skyrmion, whereas the total 
energy of the two-monopole is infinite. 


5.3 Stability of the type A, Q = 2 Skyrmopole 

A caveat of starting with a symmetric initial condition is that it may be unstable and 
only a very long relaxation time will reveal if the configuration is really stable or not. 
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(d) (e) (f) 


Figure 10: Energies for the global two-monopole inside a single Skyrmion; (a) is the 
kinetic term of the monopole (with the prefactor Q), (b) is the Skyrmion kinetic energy, 
(c) is the potential for the monopole (again with the prefactor Q), (d) is the pion-mass 
term of the Skyrmion sector, (e) is the total energy of the configuration and (f) is the 
logarithm of the total energy. The figures show slices in the xy plane for = 0. 
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Table 1: Numerically integrated energies of the type A, Q = 2 Skyrmopole. The inte¬ 
gration region is the box on which the discritized solution is found, which has the size 
203 . 


monopole kinetic energy 

16.18±0.11 

Skyrmion kinetic energy 

166.99±0.17 

monopole potential energy 

0.28±0.002 

Skyrmion potential energy 

25.49±0.03 

total energy 

208.88±0.21 


In our case of the type A two-monopole, the symmetry enjoyed by the configuration is 
axial symmetry. In order to show that the system of the two-monopole residing within 
the Skyrmion is really stable, we use again the Ansatze (5.10) and (3.3), but this time we 
perturb the axial symmetry by a finite d > 0 parameter. More specifically, we choose 


d = -ah 


(6.24) 


where a is a real number and is the lattice spacing in the x direction. We carry out 
this calculation on a 1213 with spatial-lattice spacing = hy = hz = 0.16529. 

Since the winding phases now have two centers - one at a: = —d and one at a: = d, we 
choose a profile function of the form 


h=\ tanh 


K {{x — dy + + zy 


tanh 


K ((x + dy + y^ + zy 


(5.25) 


with an appropriate value for k > 0. The fractional power of the radius in the hyperbolic 
tangent is chosen empirically to better match the d = 0 solution. In order to quantify 
the asymmetry of the configuration, we define the following spatially-weighted monopole 
charge 

Qaji = J d^x x^Q, (5.26) 

where the monopole charge density is 

Q = (5.27) 


In Fig. is shown the asymmetry between the monopole charge distribution in the x 
and the y direction for various deformation values a = 1,..., 5 as function of relaxation 
time r. In the shown relaxation-time interval, all the configurations clearly converge 
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Figure 11: Asymmetry in the monopole charge, Qx/Qy as function of relaxation time r, 
for various values of the deformation parameter o. 


towards the non-deformed configuration, i.e. d = 0. The more deformed the initial con¬ 
figuration is, the longer it takes for the relaxation to relax the configuration back to the 
non-deformed one. 


In Fig. 1^ is shown a series of snapshots of the isosurfaces of Skyrmion and monopole 
charges as the relaxation time progresses. This particular series is made with the de¬ 
formation parameter a = 5. The configuration quickly converges towards a torus shape 
and after a sulficiently long relaxation time, it becomes infinitesimally close to the a = 0 
(d = 0) solution. 

In Fig. [^is shown a series of snapshots of the contours of the norm of the monopole 
field 

+ ($ 2)2 + ($ 3 ) 2 _ ( 528 ) 

as the relaxation time progresses. This particular series is again made with the deforma¬ 
tion parameter a = 5. It is observed from the innermost contour(s) that the monopole 
field possesses two zeros when the relaxation algorithm begins. After the time of Fig. 


and before Fig. IS^, the two zeros merge to a single zero and hence becomes the non- 
deformed solution a = 0 (d = 0). We take this as a proof of stability of the torus-like 
solution shown in the last subsection. 
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; = 0,0072338 



Figure 12: The global two-monopole with a deformed initial condition (o = 5) inside 
a single Skyrmion. The figures show isosurfaces of the monopole charge and Skyrmion 
charge densities, respectively, at half their respective maximum values. The blue surface 
is the Skyrmion charge and the magenta surface is the monopole charge. As the relaxation 
time progresses, the solution converges towards that of Fig. The calculation is made 
on a 121^ cubic lattice. 
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(a) 


(b) 


(c) 




(d) (e) (f) 






(g) 


(h) 


(i) 


Figure 13: The global two-monopole with a deformed initial condition (o = 5) inside a 
single Skyrmion. The figures show contours of the norm of d>, i.e. + $2 + $3 on an 
xy slice ai z = 0. The blue circles in the middle of figures (a)-(f) show two separated zeros 
in the monopole field. As the relaxation time progresses, the two distinct zeros clearly 
merge to a single zero and the solution converges towards that of Fig. | 8 | 
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5.4 Instability of the type B, Q = 2 Skyrmopole 


Although we find that the type A, Q = 2 Skyrmopole of Sec. | 5 . 2 | is stable for the chosen 
values of the parameters in the model, we also find that the corresponding type B, Q = 2 
Skyrmopole is less stable. We carry out the same calculation as in Sec. | 5 . 2 [ but with the 


initial condition given by the Ansatze (5.10) and (3.9) with d = 0 and the parameters are 
again chosen as C 2 = 1 , C 4 = 4, m = 1 , u = 1/4, A = 128, b = 1/2 and a = 4. We take 
the profile function of the initial condition to be 


h = tanh [nr ], 

where r is the radial coordinate and k is an appropriately chosen constant. We find that 
for this choice of parameters, the type B configuration is unstable. In Fig. [T^is shown a 
series of snapshots of the isosurfaces of Skyrmion and monopole charges as the relaxation 
time progresses. We observe that the two monopole constituents split up, escape the 
Skyrmion and move off to infinity. This is, however, by no means a proof that the type B 
configuration cannot be stabilized. We suspect however, that by tuning the parameters 
to stabilize the type B configuration, it will eventually converge to the solution of type A. 


6 Conclusion 

In this paper we have shown that the two-monopole is unstable by itself both in the 
point-charge approximation and with a smooth profile function on a cubic lattice. In 
order to stabilize the two-monopole we first considered adding an interaction term to 
the monopole Lagrangian, but pointed out the difficulties in such an approach. For this 
reason we followed a different path to stabilize it. The strategy is to create an energy 
barrier for the monopoles to overcome in order for them to escape a spatial region. This 
energy barrier was created as an example by a function that multiplies the entire monopole 
Lagrangian and interpolates between a low value in the center of the region where the 
monopoles are to be contained and a high value outside. We have shown with a given 
set of parameters, the two-monopole is stable when the initial guess is of type A and by 
perturbing said configuration, we have further shown that it returns to the found solution 
after a sufficient amount of relaxation time. 

This type of stable solution has one feature which, for some purposes, may be unde¬ 
sirable. Namely, that the Skyrmion energy density dominates the monopole one on the 
scales of the charges. Of course the two-monopole has a diverging total energy as it is a 
global type of monopole, while the Skyrmion has a finite total energy. It may be possible 
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Figure 14: The global two-monopole of type B inside a single Skyrmion. The figures 
show isosurfaces of the monopole charge and Skyrmion charge densities, respectively, at 
half their respective maximum values. The blue surface is the Skyrmion charge and the 
magenta surface is the monopole charge. As the relaxation time progresses, the solution 
diverges and the monopoles eventually escape off to infinity. The calculation is made on 
a 121^ cubic lattice. 
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to find a function Q, which multiplies the monopole Lagrangian, such that the monopole 
energy locally dominates over the Skyrmion energy. The problem of the, in this paper, 
chosen type of function, is that the constraint (5.8) needs to be satisfied. A different 
function that would not have such a constraint can be constructed as 


= 2 


1 f 1 -\-b — n1 


2 + 6 


( 6 . 1 ) 


although this may be less stable than the function chosen in this paper because the energy 
barrier drops to the same value outside the Skyrmion as inside it. This type of function 
can be investigated in future works. 

One potential application of these stable multimonopole systems would be as potential 
models of dark matter halos [13 [H], as their density profiles asymptotically scale as 1/r^, 
automatically yielding fiat rotation curves. The long distance divergences also need to 
be eliminated in this context, but various mechanisms for this have been proposed in the 
literature. 


Finally, of course global multi-monopoles with charge Q > 2 may equally well be 
considered in our construction. We leave this for future investigations. 
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Appendix A The unstable two-monopole 


Eor completeness, we attempt to put the Q = 2 monopole on a lattice with the 
initial condition given by the Ansatze (5.10) and (3.3) with d = h^, i.e. the deformation 
parameter takes the values of a single lattice spacing. We take the profile function of the 
initial condition to be 

h = tanh [«;rl, 


where r is the radial coordinate and k is an appropriately chosen constant. In Eig. [T^is 
shown a series of snapshots of the isosurfaces of the monopole charge as the relaxation 
time progresses. We observe that the two monopole constituents split up, escape the 
Skyrmion and move off to infinity (at an accelerating speed). 
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(e) (f) (g) (h) 



(i) (j) (k) (1) 

Figure 15: The global two-monopole of type A with deformation parameter d = h^. 
The figures show isosurfaces of the monopole charge at half its maximum value. As the 
relaxation time progresses, the solution diverges and the monopoles eventually escape off 
to infinity. The parameters are chosen as v = 1/4 and A = 128. The calculation is made 
on a 121^ cubic lattice. 
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